Triplon mean-field theory of an antiferromagnetic model with an exact fourfold 
degenerate Shastry-Sutherland ground state 



Rakesh Kumar* and Brijesh Kumar''' 
School of Physical Sciences, Jawaharlal Nehru University, New Delhi 110067, India 

(Dated: January 21, 2013) 

In a recent work [Phys. Rev. B 77, 014419 (2008)], a quantum spin- 1/2 model having an exact 
fourfold degenerate Shastry-Sutherland ground state was constructed, and studied using finite-size 
numerical exact diagonalization. There, a schematic quantum phase diagram suggesting many 
competing phases was also proposed. In the present work, this model is investigated using the 
dimer and plaquette triplon mean-field theories. The corresponding quantum phase diagram (in the 
thermodynamic limit) is found to have many interesting phases: commensurate and incommensurate 
magnetically ordered phases, and columnar dimer, Shastry-Sutherland and plaquette spin-gapped 
phases. Among themselves, these phases undergo a continuous or a level- crossing transition. 
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I. INTRODUCTION 

The emergence of exotic physical properties in frus- 
trated quantum magnets is of great current interest.^ 
In particular, the spin-gapped systems with dimerized 
singlet ground states have been given a lot of attention. 
There are quite a few of such actively investigated ma- 
terials, of which the CuGeOa and SrCu2 (603)2 serve as 
two nice and simple examples in ID and 2D, respectively. 
Here, Cu^+ acts as a spin-1/2. Various experimental in- 
vestigations have shown that Cu-germanate is a good re- 
alization of the famous Majumdar- Ghosh (MG) quantum 
spin chain model^^, and SrCu2 (603)2 is a natural real- 
ization of the Shastry-Sutherland (SS) model. '''^ As it 
happens, the MG and the SS models are historically im- 
portant prototypical cases of how frustration and quan- 
tum mechanics lead to non-magnetic dimer (spin-singlet) 
ground states. The competing exchange interactions in 
the Majumdar-Ghosh spin chain give rise to a doubly de- 
generate dimer ground state at Ji = 2J2. The Shastry- 
Sutherland model is an extension of the MG model on 
square lattice with three-quarters of the second-neighbor 
exchange interactions deleted out. The ground state of 
the SS model is a certain exact direct product of the pair- 
wise singlets (or simply, dimers; see Fig. 2). A plaquette 
spin-gapped phase has also been shown to appear in a 
small range between the Neel and the SS phases.^' 

There are cases when the dimer ground state is not ex- 
act, but arises (or expected to arise) nevertheless because 
of strong frustration. The spin-1/2 antiferromagnetic Ji- 
J2 model on square lattice is a very good example of this. 
It is known to have two antiferromagnetically ordered 
phases: Neel and collinear. Classically, the phase transi- 
tion between these occurs at J2/J1 = 0.5. But quantum 
fluctuations reduce the sublattice magnetization, and a 
"quantum" paramagnetic region appears in the range 
0.4 < J2/J1 < 0.6.^^^^ Although there are many dif- 
ferent views on the nature of this quantum paramagnetic 
ground state, the majority opinion favors the existence of 
a columnar dimer phase. ^^'^^^^ An extended version of 
this model, the Ji — J2 — J3 model, has also been studied 



recently on square lattice around (J2 + J?) / Ji = 1/2.^^ 
A difficult question that arises in the studies of such 
problems is the nature of quantum phase transition be- 
tween the magnetically ordered (say, Neel) and the dimer 
(say, columnar) phases. For example, in the J1-J2 model, 
it is predicted by some to be a weak first order transi- 
tion. ^^^^ However, within a field-theoretic framework, 
Senthil et al proposed some years ago a generic sce- 
nario that goes beyond the Landau- Ginzburg- Wilson 
paradigm. ^^'^^ They proposed the possibility of a con- 
tinuous quantum phase transition between a dimer and 
a magnetically ordered phase without using order param- 
eter notions. Instead, the concept of "deconfined" quan- 
tum criticality in which free spinous at the quantum crit- 
ical point act as the essential degrees of freedom is put 
forward. Motivated by this remarkable suggestion, many 
quantum spin models have been constructed and inves- 
tigated. ^^'^^'^^ The signatures of a deconfined quantum 
critical transition seem to have been observed in a SU(2)- 
invariant spin-1/2 multiple-spin exchange model on the 
square lattice. However there is no finality yet. There 
are, for example, counter suggestions wherein a ring- 
exchange model, which breaks the Landau-Ginzburg- 
Wilson framework, shows a first order phase transition. 

With similar motivations, Gelle et al. constructed and 
studied a quantum spin-1/2 model on square lattice a few 
years ago.'^^ The important part of this model is a certain 
projector interaction term (leading to four-spin exchange 
interactions; described in the following section) which 
helps in realizing the SS dimer state as exact ground 
state. Interestingly, this SS ground state is fourfold de- 
generate (unlike in the original SS model), as it sponta- 
neously breaks the lattice symmetry. Adding the nearest 
or further neighbor Heisenberg interactions to this ex- 
actly solvable model generates a competition to the ex- 
act SS ground states, and making it more interesting as 
well as difficult to study. They investigated it by ex- 
act diagonalization and presented a schematic quantum 
phase diagram. In this paper, we study this model using 
triplon mean-field theory with respect to various dimer 
and plaquette phases, and present a rich quantum phase 
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diagram. 




FIG. 1. (Color online) This figure represents pictorial rep- 
resentation of the Hamiltonian Eq. (1). Filled circles and 
squares denote the spins on two sub-lattices. Broken dotted 
lines represent Ji interactions whereas solid lines J2 inter- 
actions. Similarly, J3 interactions can be understood (not 
shown here) . The shaded triangles with opposite orientations 
represent Hk- 



II. MODEL 



The operator P3/2 annihilates those states for which the 
total spin of the triangle is 1/2. Otherwise, it is one. 
For this particular choice of the projector interaction, 
the Hamiltonian (in the absence of the Heisenberg in- 
teractions) was shown to have an exact spontaneously 
dimerized singlet ground state, consisting of four degen- 
erate Shastry-Sutherland configurations (see Fig. 2). For 
the details of the exact solution, see Ref. 31. 
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We consider the following spin- 1/2 model on the square 
lattice. 



H — Ji ^ • Sj + J2 ^ • Sj 



J3 



E 



{iJ) 



+K 



E 

[i,j,k;l,'m,n] 



{{iJ)) {{{id))) 

T^3/2{iJ^k)V3/2{l,rn,n). 



(1) 



The couplings Ji, J2, and J3 are respectively the first, 
second, and third neighbor antiferromagnetic Heisenberg 
exchange interactions. Originally in Ref. 31, only Ji 
and J2 were considered. But we have also included J3 
for the reasons that will become clear shortly. The last 
term of strength K is the projector interaction. Here, 
K > 0, and the summation [i, k;l,m,n] denotes the 
sum over all horizontal and vertical (six-spin) plaquettes 
(see Fig. 1). On each plaquette, we designate two oppo- 
sitely oriented triangles with vertices and /,m,n. 
The projector interaction is thus given by the product of 
the spin-projection operators, 7^3/2, defined on the two 
triangles of a plaquette. The projection operator 7^3/2 
on a triangle is defined below. 



7^3/2 ihj^k) 



(S,-S, 



Sj • Sk 



Sk-Si) (2) 



FIG. 2. (Color online) The four SS dimer states. These states 
spontaneously break the lattice symmetry and connected by 
lattice translations in x and y directions. 

The projector interaction, when expanded, would have 
the two-spin (Heisenberg) terms, and certain four-spin 
exchange interactions. Since the two triangles of a (ver- 
tical or horizontal) plaquette belong to two different 
sublattices of the square lattice, only J2 and J3 type 
Heisenberg exchange would arise out of the K interac- 
tion. Therefore, it seems natural to also consider the 
competition from J3, in addition to Ji and J2 as done in 
the original model. 

Using the definition of projection operator in Eq. (2), 
the Hamiltonian Eq. (1) can be rewritten as 

n = ^KL + Ji ^ s, • s, 

+J2 Si • Sj + J3 Si • Sj 

{{iJ)) {{{iJ))) 

^-K ^ h{i,j,k)h{l,m,n), 

[i,jf,/c,Z,m,n] 

where L is the number of lattice sites, J2 
-^3 = -^3 + I and /i(z, j, k) = Si- S 



(3) 
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Sj • Sk ■ 



3 ' 
- Sh • Si. 
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In the following sections, we will discuss this model using 
triplon mean-field theory. 



1. On Shastry- Sutherland dimer state 



III. TRIPLON MEAN-FIELD THEORY 

As shown in previous section, the Hk Hamiltonian 
gives four symmetry broken SS dimer states. This means 
that the model is exactly solvable at the limit Ji = J2 = 
J3 = 0. The triplon mean- field theory {also known as 
bond- operator mean- field theory) is a very good analytical 
approach on the situations where spontaneous dimeriza- 
tion of singlet dimers happen as in the present case. We 
formulate triplon mean- field theory on a chosen SS dimer 
state. Moreover, The previous research of J1-J2 model on 
square lattice predict columnar-dimer^^'^^'^^, plaquette- 
RVB^^^^ states in the paramagnetic region between Neel 
and collinear phases. So, we also devised the triplon 
mean-field theory on columnar dimer and plaquette-RVB 
states. 

In this section, we describe two kind of triplon mean- 
field theories: (1) Dimer triplon mean- field theory for SS 
and columnar states, and (2) Plaquette triplon mean-field 
theory for plaquette-RVB state. 



A. Dimer triplon mean-field theory 

The Hilbert space of a spin dimer is four dimensional 
that is one singlet \s) and three triplet \tx)^ \ty)^ \tz) 
states. We transform these states into Fock space. In 
the Fock space, the singlet and triplet states are forming 
by the application of creation operators on the vacuum 
|0). These operators are known as bond-operators^^ and 
defined as 
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By calculating matrix elements of Si and S2, we can 
express the spin operators in terms of bond operators as 



5f 
5? 



-tU 



(5a) 
(5b) 



where e^/Sj is the totally antisymmetric tensor. The bond 
operators are canonical bosons and satisfy commutation 
algebra. The number of bosons on a dimer is restricted 
by the constraint: s^s + ta^a = 1, where a = x^y^ z. 




FIG. 3. (Color online) Figure (a) shows the Shastry- 
Sutherland dimer state, figure (b) is its unit cell. The cor- 
responding Brillouin zone is shown in figure (c). 

We choose a reference SS dimer state shown in Fig. 3. 
For the sake of generality, we assume different singlet 
condensation on each dimer in unit cell. The singlet con- 
densation amplitudes on the dimers (1) and (2) are some 
c-numbers Si and 52 respectively. We also used distinct 
type of triplet operators on these dimers, say ^ and 

(2) 

the ' The constraints are imposed by introducing global 
chemical potentials /i^^^ and /i^^^ on the dimers of unit 
ceh: 




+M2 



.-^+^t«t(R)tW(R)-l 

a 

s-^+^e)t(R)e)(R)- 



(6) 



where R refers to position of unit cell. The chemical 
potentials fix the number of bosons on the dimers. We 
ignored the triplet-triplet interactions as they change the 
results very minutely. Using these approximations, the 
mean-filed Hamiltonian can be written as 



(7) 



where. 



^0 



rrk 



-L 



-J2(5l + s1)+A«S? + A(2)5i" 

k,a i=l,2 



(8) 
(9) 



^-ka^-k 



The renormalized effective chemical potentials A^*^ and 
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^^"^ are give by following equations 



A(0 



J2 - /i 



(10) 



J2 (cos + cos 2ky + cos 2kx cos 



— J3(cos + cos 2ky) — -Ksl cos 2kx cos 2A:^11) 

o 



A close inspection of Eqs. (9) tells us that the mean- 
field Hamiltonian decoupled into independent triplet op- 
erators that is no mixing of triplet operators t^^ and 
t'a'- So, we can further simplify Eqs. (8), (9), (10), 



and (11) by taking si = S2 



(1) 



(1) 



-5, i'a 



t 



(2) 



try, and 



/i. We diagonalize the resultant mean- 
field Hamiltonian by unitary Bogoliubov transformation: 

t]^a = 7ka cosh6>k + T^ka sinh6>k, where 7ka are Bogoli- 
ubov bosons. The diagonalized mean-field Hamiltonian 
thus obtain is following 



U 



mf 



(12) 



In the above equation, = ^^A (A — 2s^^k) > is 
triplon quasiparticle dispersion. The energy of the mean- 
field Hamiltonian Eq. (12) gives ground energy in the ab- 
sence of quasi particles. The saddle point equations of the 
ground state energy with respect to unknown mean-field 
parameters 5, A lead to the self-consistent equations. The 
solution of these equations give mean-field results. The 
detailed procedure of finding solution of self-consistent 
equations have been discussed in Ref. 19. 



2. On columnar dimer state 

In the case of J1-J2 model on square lattice, the triplon 
mean-field theory predicts columnar dimer state in the 
range 0.19 < J2/J1 < 0.67 (see Ref.^'^). The present 
Hamiltonian also has these two interactions, Ji and J2. 
In this regard, we formulate triplon mean-field theory on 
a reference columnar dimer state shown in Fig. 4. This 
dimer state has single dimer per unit cell. We assume 
singlet Bose condensation on the dimers (shown thick 
line segments in Fig. 4). Let say it is some c- number s. 
Here again, we ignore triplon interactions and the con- 
straint, + Xla^a^Q! = is applied globally by choos- 
ing chemical potential ja on each bonded spins uniformly: 
[^^ + Z]q; ^a(I^)^Q;(I^) ~ 1] 5 where R denotes posi- 
tion of a dimer. The triplon mean-field Hamiltonian of 
the model on the CDS state is now can written as 



U 



mf 



Eq 



rrk 



(13) 




ky 

(7r/2,7r) 



(a) 



(b) 



FIG. 4. (Color online) The figure (a) is columnar dimer state. 
The thick line segments are singlet dimers. Figure (b) show 
the Brillouin zone of the columnar dimer state. 



where 



k, a 



-ko! 



(14) 



(15) 



In the above equations, A = | Ji — /i is effective chemical 
potential and ^/^ is given by 

= -^i(cos 2kx — 2 cos ky) + J2(cos 2kx cos ky + cos ky) 



- J3 (cos 2k X -\- cos 2ky). 



(16) 



The mean-field Hamiltonian is diagonalized by unitary 
Bogoliubov canonical transformation. The diagonalized 
mean-field Hamiltonian can now be written as 



k, a 



(17) 



The dispersion relation is given by = \/ A(A — 2s^^]^. 
Here, 7 are Bogoliubov quasi-excitations. 



B. Plaquette triplon mean-field theory 

Consider the labelled plaquette shown in Fig. 5. The 
antiferromagnetic Heisenberg Hamiltonian on plaquette 
can be written as 



Hp — Ji(Si • S2 + S2 • S3 + S3 • S4 + S4 • Si) 

+J2(Si- S3 + 82-84). 



(18) 



Let the total spin on bonds (13) and (24) are S13 and 
S24 respectively. Then, the total spin of the plaquette 
will be 

S = Si3 + S24; (Si3 = Si+S3 and S24 = S2 + S4). 

(19) 
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FIG. 5. (Color online) The plaquette resonating valence bond 
state. The labelled plaquette is chosen as the unit in the 
plaquette triplon mean-field formulation. 



TABLE I. The eigen spectrum of the plaquette Hamilto- 
nian 1-Lp. In the table, the singlet and triplet [Sz — 0) 
state are denoted by [i^j] and {i^j} which are defined as 

[hj] = 75 (I t4.) + 1 u^j)). {hj} = ^^{\ t4.> - 1 ;.t.». 

5'i3 6*24 S Eigen energy 




1 

2 



6 J a 
"2 ^2 

"2 

-2Ji + |J2^ 
— Jl + 2^2^ 
Jl + 1^2^ 



^ Eigenstate: [1, 2] (g) [3, 4] - [2, 3] (g) [4, 1]. 
^ Eigenstates: 

[1, 3] ^ I UU), [1, 3] (8) ^ (I ^2U) + I ;2t4» , [1, 3] (g) I t2t4>. 

^ Eigenstates: 

I ;i;3> ® [2,4], ^ (I ti^s) + I ;it3» ® [2,4], I tit3> ® [2,4]. 
^ Eigenstate: ^ ([1, 2] (g) [4, 3] + [1, 4] (g) [2, 3]). 
^ Eigenstates: 

^ I ^i^s) ® {2, 4} + ^{1, 3} ^ I ;2;4>, ^ I tits) ^ I ;2;4>+ 

I ;i;3> ^ I t2t4>, -^1 tlt3> ® {2, 4} + ^{1, 3} ^ I t2t4>. 

^ Eigenstates: 

I ;i;3> (8) I ;2;4>, ^1 uu) ^ {2,4} + ^{1,3} (g) I ;2;4>, 
^{1, 3} ^ {2, 4} + I ;i;3> ^ I t2t4> + ^ I tits) ^ I ;2;4>, 

^\ tlt3> ^ {2,4} + ^{1,3} ^ I t2t4>, I tits) ^ I t2t4> . 



Using Eq. (19), we can rewrite Eq. (18) as follows 
Hp = ^(S^ - - S^4) + J2(Si • S3 + S2 • S4) . (20) 



The obtained energy spectrum of the Hamiltonian Hp 
is given in Table (I). If Ji is greater than J2, we can 
restrict ourselves in the restricted subspace of the singlet 



state with energy — 2Ji + | J2 and the three degenerate 
triplet states having energy — Ji + ^ J2 . These states 
can viewed as the states emerging out of vacuum by the 
boson creation operators as \s) := s^|0), \t+) := t\_\0), 
\t^) := 4|0) and := 1^0), where 

V ([1,2]® [4, 3] + [1,4]® [2, 3]), (21a) 



\t-) 



V3 

^1 tits) ® {2, 4} + 3} ® I t2t4), (21b) 

^1 tits) ® I t2t4) + ^1 tits) ® I t2t4), (21c) 

^1 tits) ® {2, 4} + 3} ® I t2t4). (21d) 



The matrix elements of the spin-operators, Sx^ Sy^ and 
Sz, in the subspace 1^^), give the following 

definition of the spin-operators 

i-Y' 



(4s + sH,) + ^ (tlt+ - tU.) , (22a) 
= ^ (^+« + «^*-) + Ifjn (^+*- - 4*-) , (22b) 



\/3 



^ (tls + sH+) + ^ (t\t+ - tlu) , (22c) 

where the vertices of the plaquette are labelled by m. To 
fix the boson numbers on each plaquette, the constraint 
s'^s + t^j^t+ + t\tz + t^-t- = 1 is applied therein. Let 
introduce the operators tx and ty as following 

tx = ^{t++t-), (23) 



V2 



(24) 



Now, the constraint equation becomes s^t\.tx -\-t^yty + 
tltz = 1, and the spin-operators shown in Eq. (22) can 
be written as 



(-r 

\/6 



(25) 



where o^, /3, 7 = x^y^z and eQ,^^ is Levi-Civita antisym- 
metric tensor. 

In the simplified triplon mean-field theory, we ignore 
the mixing of different triplet bosons and higher order 
triplet interactions. Then, the exchange interactions 
within plaquette and inter-plaquette can be written as 
following 

Up (R) = f-2Ji + 5t(R)s(R) 



+ f-^i + J^2)Ett,(R)t«(R), (26) 



2 



^5;^(R)5;?(R') « LJ_,-2^ [i^(R)4(R') 
+4(R)t«(R') + t„(R)t«(R') + 4 (R)tt (R')] ^ (27) 
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where R and are the position vectors of different pla- 
quettes. 

Let the singlet Bose condensation amphtude on 
the plaquettes (see the shaded plaquettes in Fig. 5) 
is some c- number 5, i.e. (s^) = (s) = s. 
The global constraint is imposed on each plaque- 
tte: /i [s^ + tl(R)ta{'R) - 1] . The mean-field 
Hamiltonian now becomes 

Hmf = Eo ^ ^ I^Ak (^^ka^kQ, + ^ka^ka 



k, a 



+Ak ^^ko!^ — ko; H~ ^kcn^ — kcK 



(28) 



where 
L 



Ak 
Ak 



-K 



5x 2 



-K 



12 



1 

36^ 



-Ks' 



- Ji + J2 + 2 J3 



-Ks' 



(cos 2kx + cos 2/c^) + J2 cos 2A:a; cos 2/c^ 



(29) 

(30) 

(31) 
(32) 



In the above equations, we have used the defini- 
tion of renormalized effective potential, that is, A = 

^ (y—Ji + ^^2 — We diagonalize the mean-field 

Hamiltonian Eq. (28) by Bogoliubov transformation. 
Thus, we get the following diagonalized form of the mean- 
field Hamiltonian 



(33) 



where E\^ = y^^A^^^A^ > is the triplon dispersion. 



IV. RESULTS AND DISCUSSION 

In this section, we discuss the results for three different 
cases of the model. In these cases, parametrization of the 
exchange couplings is as follows 

Case (a): Ji — J2 — K model 

Ji = (1 - C) J, J2 = CJ,J3 = 0,K = i^,J^f^ = l, (34) 
Case (h): J2 — J3 — K model 

Ji = 0, J2 = (1 - 7?) J, J^ = r]J,K = ^- n = l, (35) 
Case (c): Ji — Js — K model 

Ji = (1 - 5) J, J2 = 0, J3 = ^ J, = J + = 1, (36) 

where all introduced parameters are positive as the cou- 
plings Ji, J2, J3, and K are antiferromagnetic. In the 
following discussion, we focus on these three cases. 

Case (a): In this case, we extract physics from dimer 
triplon mean-field theory as follows. First, we solve the 



self-consistent equations obtained by dimer triplon mean- 
field theory on SS dimer state. It gives the phase bound- 
ary between collinear phase and SS dimer state. Im- 
portantly, the SS dimer state does not compete with 
Neel state as the SS dimer state is also an eigenstate for 
the nearest-neighbor antiferromagnetic Heisenberg inter- 
action, Ji -s^ Si • Sj . Thus, we get an overestimated 
phase boundary of the collinear state. We can evade from 
this overestimation by formulating the triplon mean-field 
theory on the columnar dimer state (CDS). This is phys- 
ically correct adoption as, for Ji — J2 model, the triplon 
mean-field results show columnar singlet dimerization in 
the range 0.19 < J2/J1 < 0.67 The triplon mean- 

field analyses on the CDS show two commensurate mag- 
netically ordered phases: Neel and collinear, and an in- 
commensurate phase with ordering wave vector (0,(/?), 

where cp = cos~^ (2^)* triplon mean-field the- 

ory on the plaquette crystal state gives an energetically 
more favorable PRVB phase inside the columnar dimer 
phase. The resultant quantum phase diagram is shown 
in Fig. 6. We explored plaquette valence-bond crystal 
phase in the quantum phase diagram because some ear- 
lier studies show the existence of plaquette phase in the 

paramagnetic region of the Ji — J2 model on square lat- 
tice^o 22, 




FIG. 6. (Color online) The quantum phase diagram for the 
J1-J2-K model. The solid black lines denote continuous phase 
transitions whereas broken lines and small circles indicate the 
first-order transitions. In the phase digram, we have three va- 
lence bond crystal phases: CDS, SS, and PRVB. The ordering 
wave vectors of the Neel and collinear phases are respectively 
(0,7r) and (0,0). The incommensurate ordered phase has or- 



dering wave vector (0, (p), where = cos ^ 
ing wave vectors are with respect to CDS. 



All order- 



The phase diagram contains continuous as well as first- 
order phase transitions. The Neel phase shows contin- 
uos phase transition with columnar dimer phase. The 
collinear phase also has continuous quantum phase tran- 
sitions with columnar and SS dimer states. The incom- 
mensurate magnetic phase shows first order phase tran- 
sition (i.e., level crossing) with SS dimer phase and con- 
tinuos phase transition with columnar dimer phase. The 
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paramagnetic phases, CDS and SS, show first order quan- 
tum phase transition with each other. The plaquette 
valence bond crystal phase shows a level crossing with 
CDS. The exact diagonalization study'^^ shows the phase 
boundary of the collinear phase upto = 0.87 at J2 = 1 
axis whereas our study underestimating and giving 0.52. 
At Ji = 1 axis, the Neel phase is bounded by < 0.2 
(exact diagonalization), 0.08 (triplon mean- field theory). 
Again, triplon mean-field theory underestimating. How- 
ever, the SS dimer phase exists for > 0.58 at Ji = 1 
axis in triplon mean-field theory which is very close to 
numerical value, 0.6, obtained by exact diagonalization 
on 32 site cluster. The numerical study predicts many 
kind of phases in between the Neel phase and the SS 
dimer phase. This is happing because of the presence 
of different competing dimer correlations as revealed by 
triplon mean-field analyses. 

When we are close to a magnetically ordered phase, 
the triplon mean-field results vary from exact numerical 
results. For example, in the case of Ji — J2 model, the 
exact diagonalization predicts the paramagnetic region 
in the range 0-4 < J2/J1 < 0.6 11^^ whereas the triplon 
mean-field theory gives 0.19 < J2/ Ji ^ 0.67 i'^'^^. This is 
because we always assume singlet condensation in triplon 
mean-field theory calculation and triplons are created on 
the top of the frozen singlets. This means that there is 
always phase coexistence in ordered phase, that is, {s) 7^ 
and {to,) / 0. 
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FIG. 7. (Color online) The quantum phase diagram for the 
J2-J3-K model. The ordering wave- vectors on the SS dimer 
unit cell. In this phase diagram, there are two magnetically 
ordered phases: collinear phase and an incommensurate phase 
with ordering wave vectors (0,0) and (7r/2,7r/2) respectively. 
The ordering wave vectors are with respect to SS dimer state. 

Case (b): In the absence of Ji interaction, the model 
consists only second and third nearest-neighbor interac- 
tions. This leads to competition between a antiferromag- 
netic phase and SS dimer phase. The quantum phase di- 
agram for this case shown in Fig. 7. This is obtained by 
formulating the triplon mean- field theory on the SS dimer 
state. Here, we have two magnetically ordered phase with 
ordering wave vectors (7r/2,7r/2) and (0,0) (i.e., collinear 




K 



FIG. 8. (Color online) The quantum phase diagram for the 
J1-J3-K model. The solid black lines denote the continu- 
ous transitions and broken line represents first-order quantum 
phase transition. The (0,7r) phase is Neel phase. The incom- 
mensurate phase has ordering wave vector (7r/2, tt — (/?), where 

cp = cos~^ (4^)" wave- vectors defined with respect 

to CDS. 

phase). The SS dimer state gives continuous quantum 
phase transitions with collinear and (7r/2,7r/2) magneti- 
cally ordered phases. 

Case (c): For this case, the quantum phase diagram is 
shown in Fig. 8. We do the triplon mean- field analyses on 
the columnar and SS dimer states. The phase diagram 
contains Neel phase, an incommensurate phase with re- 
spective ordering wave vectors (0,7r) and (7r/2,7r — (/?), 

where = cos~^ (4^)* dimer state shows en- 

ergy level crossings with the incommensurate phase and 
columnar dimer phase. The columnar dimer phase shows 
continuous quantum phase transitions with Neel and in- 
commensurate phase. 

V. SUMMARY 

We studied the model proposed by Gelle et al.^^ using 
triplon mean-field theory with respect to various dimer 
and plaquette phases. We presented the quantum phase 
diagrams for the three different cases of the extended 
model: (a) Ji — J2 — K model (6) J2 — J3 — K model, and 
(c) Ji — Js — K model. In our analyses, we used dimer 
and plaquette triplon mean- field theories. In Ji — J2 — K 
model, we get three valence bond crystal phases namely 
SS, CDS, and PRVB. There are three magnetic ordered 
phases emerge out with ordering wave vectors (0,7r), 

(0,0), and (0, (/:?), where (p = cos~^ (2^)' phase 
transition among these magnetic ordered and valence 
bond crystal phases are either continuous or first-order. 
In the case of J2 — J3 — K model, the commensurate 
ordered phases (0,0) and (7r,7r) show continuous phase 
transitions with SS dimer state. In the last case, that 
is for Ji — Js — K model, the CDS phase shows contin- 
uos phase transition with Neel and an incommensurate 
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(7r/2,7r - if) phases, where if = cos"^ (4^)- 

dimer state shows first-order transitions with CDS and 

ordered incommensurate phases. 
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